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Abstract 

We present a simple model which describes elastic response of single-stranded DNA 
(ssDNA) to stretching, including the regime of very high force (up to 1000 pN). 
ssDNA is modelled as a discreet persistent chain, whose ground state is a zigzag 
rather than a straight line configuration. This mimics the underlying molecular 
architecture and helps to explain the experimentally observed staturation of the 
stretching curve at very high force. 
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1 Introduction 



Over the recent years, impressive progress has been made in understanding 
of the physics of DNA. To a large extent, this development was due to the 
introduction of novel micromanipulation techniques suitable for study of the 
mechanical response of a single molecule [1J[2J. The mechanical properties of 
double-stranded DNA (dsDNA) have been studied in pioneering experiments 
of Smith et al pQ. These experiments agree remarkably well with Worm-Like 
Chain (WLC) model which describes DNA as an elastic rod, whose bending 
modulus is proportional to chain persistence length l p (l p w 50nm for dsDNA, 
which is the only adjustable parameter of the model) [3]. 

In the case of single-stranded DNA (ssDNA) and RNA molecules the choice 
of an adequate elastic description remains an open problem. On the one hand, 
ssDNA and RNA were traditionally described by Freely Joint Chain model 
(FJC). Its extensible version has been originally used for fitting of the early 
ssDNA stretching data [2]. On the other hand, the chemical structure of ssDNA 
strongly support the picture with a finite bending modulus, reminiscent of 
WLC model. However, being a continuous model, the WLC description is 
unlikely to be valid in the regimes when discrete nature of chemical bonds 
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becomes relevant (e.g. for a sufficiently high stretching force). To overcome 
these limitations, Discreet Persistent Chian (DPC) model has been proposed 
for ssDNA by Storm and Nelson [3]. Interestingly, WLC itself was introduced 
as a continuous limit of a similar discrete model, originally proposed back 
in late 40s by Kratky and Porod (KP) [5J. Authors of Ref. [I] succeeded in 
calculating the response of DPC to an arbitrary stretching force, by combining 
the Transfer Matrix method with a variational procedure. Nevertheless, the 
discreet model yields only a marginal improvement of the fit of ssDNA force- 
extension curve, compared to WLC and FJC cases. 

In order to improve the fitting, the extensible versions of the above models 
are typically used, i.e. the linear deformations of the bonds are introduced. 
However, the AFM experiments done in the limit of very high force (up to 1000 
pn) [B] , exhibit surprising saturation of the end-to-end distance at certain value 
above the chain counter length. Any model which treat deformations of the 
ssDNA backbone in a linear manner, is unable to capture this phenomenon. 
The experimental behavior is however consistent with results of detailed ab- 
initio calculations [7]. In this paper, we propose a simple model motivated 
by actual microscopic architecture of ssDNA. Without introduction of extra 
fitting parameters, our "zigzag" model adequately describe both low- and 
high-force regimes of ssDNA stretching. 



2 Zigzag model for ssDNA elasticity 



Discreet Persistent Chian (DPC), also known as Kratky and Porod (KP) 
model, has the following Hamiltonian: 



J -®l l+l - bf ■ t, 



Here tj, is i-th. bond orientation, b is (fixed) bond length, 6^+1 is the angle 
between i — th and i + 1-th bonds (i.e cos 6^+1 = tjtj + i), and f is external 
stretching force. Interestingly, two other common models can be obtained 
from DPC Hamiltonian, by taking limits of zero bond length (WLC), or zero 
coupling J (FJC). 

At the atomistic level, the backbone conformation of ssDNA is often parame- 
terized by several torsional angles per nucleotide, as shown in Figured! From 
this point of view, WLC and DPC models describe small deviations from the 
ground state configuration. . However, at the microscopical level, all the valent 
angles have certain preferred values different from 180°. The minimal general- 
ization of DPC model which mimics this observation is a chain whose ground 
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Fig. 1. Microscopic architecture of ssDNA. 
state is zigzag (all-trans) configuration: 



H = 2 + ©o cos T/>i,» + i) - 6F • tj 



Here ipi is the torsional angle associated with rotation of i — th bond. For sim- 
plicity, we assign the same preferred angle O to each bond pair. Surprisingly, 
this minimal model results in a significant improvement over the previous at- 
tempts to describe the experimentally observed large force behavior of ssDNA. 



3 Mapping on DPC model 

In order to calculate the stretching curve of the "zigzag" chain, we note that 
the model can be mapped back onto DPC model, with renormalized param- 
eters. Let us consider a chain consisting only of odd (or only even) nodes of 
the zigzag (i.e. "i" — "2 + 2" — "i + 4" —...). The ground state of this chain 
is a straight line, but the bond length (even without thermal fluctuation) is 
a function of applied force: b* = 2bcos6*, where 9* (f) is the preferred bond 
orientation angle 9i, with respect to the direction of f. In ground state con- 
figuration (i.e. zigzag) 9- t = (— 1)* 9* The value of 9* may be obtained by 
minimizing the Hamiltonian, at given stretching force: 



H ( 9 *) = E 2 (2r " 00)2 ~ fbcos6 * 



(3) 



In the limit of small Gq, we obtain:: 
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It is also easy to demonstrate that the renormalized coupling in the final 
DPC Hamiltonian will be J* = J/2. The parameters of the renormalized 
model are related to the experimentally observable persistence length as l p = 
b* (0) J*/kT fh Jb/kT. We can now use DPC model to calculate the stretching 
curve z/L = x (/), and include the effect of renormalization in the following 
manner: 



z 
L 



b* (f) 
b* (0) 



*(/) 



( 



\ 



1 + 



I, 



f}y_ 

2kT 



B 2 



*(/) (6) 



Here z and L are end-to-end distance and the original chain counter length, 
respectively. Since the discreetness of DPC is only marginally important, one 
can relate x to / with a simple interpolation formula proposed by Marko and 
Siggia for continuous WLC model: 
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Fig. 2. Stretching behavior of ssDNA, given by the zigzag model. Experimental 
points are from Refs [6] (circles) and [2] (diamonds) . Insert: Fitting of the experi- 
mental data with extensible FJC, WLC and DPC models [3]. 

The resulting stretching curve is shown in Figure [2j It is remarkable that the 
very good fit can be achieved with the same number of free parameters, as in 
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the earlier models: 6(which is fixed at its physical value of a typical chemincal 
bond length, O.lnm ), ssDNA persistence length l p 0.85nm (extracted from 
the moderate force stretching behavior). The only fitting in our model is angle 
6 /2 fa 25° which essentially replaces the bond rigidity commonly used in 
the extensible versions of classical models. 



4 Conclusions 

In conclusion, we proposed a new model for ssDNA elasticity. It differs from 
the earlier ones by an explicit account for certain features of microscopic ar- 
chitecture of the chain backbone. Qualitatively, our model attribute the ultra- 
high-force stretching behavior to the perturbations to the valent angles rather 
than deformations of the bonds themselves. This implication can be directly 
checked by the first principle calculations. 
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